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ABSTRACT 

Dominated by a single, large- amplitude pulsation mode, the rapidly-pulsating hot 
subdwarf B star CS 1246 is a prime candidate for a long-term O-C diagram study. 
We collected nearly 400 hours of photometry with the PROMPT telescopes over a 
time span of 14 months to begin looking for secular variations in the pulse timings. 
Interestingly, the O-C diagram is dominated by a strong sinusoidal pattern with a 
period of 14.1 days and an amplitude of 10.7 light-seconds. Underneath this sine 
wave is a secular trend implying a decrease in the 371. 7-s pulsational period of 
P = —1.9 x 10~ n , which we attribute to the evolution of the star through the H-R 
diagram. The sinusoidal variation could be produced by the presence of a low-mass 
companion with msini ~ 0.12Mq orbiting the subdwarf B star at a distance of 20 
Rq. Analysis of the combined light curve reveals the presence of a low-amplitude 
first harmonic to the main pulsation mode. 

Key words: stars: individual: CS 1246 - stars: oscillations - stars: subdwarfs 



1 HOT SUBDWARFS AS BINARY SYSTEMS 

Subluminous B (sdB) stars represent one of the least- 
understood stages of stellar evolution, yet they dominate 
surveys of faint blue objects and are found in almost all 
Galactic stellar populations. They are the field counter- 
parts of the Extreme Horizontal Branch (EHB) stars in 
globular clusters and lie in a region of the H-R diagram 
corresponding to the location of model s with He-bu rning 
cores and extremely thin H envelopes (|Heberlll986l ). Un- 
derstanding these stars better would illuminate the enig- 
matic "second parameter" problem in globular cluster evo- 
lution, improve synthesis modeling of the UV upturn in 
giant elliptical galaxies, identify the evolutionary mecha- 
nisms responsible for guiding main sequence stars to this 
state, and constrain the physical properties of the dense 
plasmas present in these compact objects. 

Presumably, the sdB stars have been stripped of al- 
most all their surface hydrogen after the red giant phase, 
but how and why this stripping occurs is still debated. 
Some have proposed binary models in which the extra an- 
gular momentum provided by a companio n manages to ex- 
pel a lmost all of the hydrogen envelope l|Han et al.l [2002. 
l200Sf ) . Others have proposed that more modest amounts of 



angular momentum, as migh t be resident in a planet, could 
accomplish the same effect (|Sokerlll998l l. Still, others are 
able to draw hot subdwarf stars from their models w ithout 
the assistance of a companion ijD'Cruz et al.l[l996u . 

In any case, companions may be profoundly important 
for understanding the evolution of sdB stars, and placing 
constraints on the nature of these systems may shed light 
on their evolutionary histories. Several studies have shown 
an unusuall y high binary fracti o n amongst the hot sub - 
dwarf stars dMaxted et al.ll200ll ; iNapiwotzki et al.ll2004 ). 
Measurements of radial velocity (RV) shifts have led to the 
orbital paramet ers of approximate ly 85 sdB binaries (see 
Table A.l from iGeier et al.ll2011dY The periods of such 
systems range from 0.07 to 28 days, with a peak between 
0.5 d & 1.0 d. RV variability studies are most sensitive to 
systems with shorter periods and higher total masses, and 
as a result, few system parameters exist for longer-period 
binaries. Only two systems with periods above 10 days 
have well-determ ined orbital parameters (P G 0850+170 
& PG 1619+522: [Morales-Rueda et al.ll2003h . 

T he discovery of the firs t of the pulsating sdB (sdBVQ) 
stars dKilkennv et al. | ll997l) in conduction with the first 
sdBV star models ijCharpinet et al.1 Il996l . Il997h opened 
the possibility of exploring their structure through a com- 



* Based on observations at the SOAR Telescope, a collabora- 
tion between CPNq-Brazil, NOAO, UNC, and MSU, and on 
observations collected at the PROMPT array via SKYNET. 
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servatories, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc, under contract with the 
National Science Foundation 



1 Her e we use the nomenclature proposed by iKilkennv et alj 
fcOlOh . All sdB pulsators are given the basename 'sdBV; sub- 
scripts V and V are added upon observations of rapid and 
slow pulsations, respectively. For reference, the sdBVv stars 
have also been called V361 Hya and EC 14026 stars, while the 
sdBV s stars have the aliases PG 1716, V1093 Her, and 'Betsy' 
stars. 
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bination of spectral analysis and stellar seismology. Multi- 
colour photometry and time-series spectroscopy tech- 
niques, for example, have revealed information about the 
total mass, envelope mass, and chem i cal st ratification for 
a handful of pulsators (see l0stensenl ((2009) for a review). 
The pulsations may also provide the opportunity to detect 
companions by using phase measurements t o look for os- 
cillati ons in the arrival times of the pulses. ISilvotti et al.l 
(|2007l ) used this method successfully to detect a 3.2Mj up 
planet around the star V391 Peg (HS 2201+2610). Their 
study also resulted in the first published measurement of 
P for a pulsating sdB star, constraining the evolutionary 
state of V391 Peg. 

CS 1246 , a recently-discovered sdBV r star 

(|Barlow et all 120101 ). is a prime target for an ex- 
tended O-C diagram study since it exhibits a single, 
large-amplitude pulsation mode. Using the PROMPT 
telescope array, we monitored the pulse timings for 14 
months by collecting nearly 400 hours of time-series pho- 
tometry. The O-C diagram shows a sinusoidal oscillation 
superimposed on a parabola. We interpret the quadratic 
term as a secular decrease in the pulsation period, 
while the sinusoidal oscillation implies the presence of 
a low-mass stellar companion, probably an M-dwarf or 
white dwarf. Our study marks the second time P and 
the mass of a companion have been estimated for a hot 
subdwarf star using the O-C method. 



2 PHOTOMETRIC DATA 

We began monitoring CS 1246 in 2009 using the 0.41-m 
Panchromatic Robotic Optical Monitoring and Polarime- 
try Telescopes (PROMPT) on Cerro Tololo to look for 
phase variations in the pulsations. The PROMPT operate 
under the control of a prioritised, queue-scheduling system 
called SKYNET and are 100% automated. Priority levels 
are assigned to all observations in the queue, and, conse- 
quently, some of our photometry runs were interrupted by 
higher-priority targets such as gamma-ray bursts. For ad- 
ditional detai ls on the PROMPT a nd SKYNET, we refer 
the reader to lReichart et all (|2005T ). 

Over the course of 14 months, we collected more than 
33,000 individual frames representing 400 hours of time- 
series photometry. The majority of these data were taken 
in a four-month span between January and May 2010 with 
PROMPT 3. On a typical night, we obtained at least three 
hours of photometry with 30-s exposures and no filter; dur- 
ing these observations, our duty cycle was on average 83%. 
We also incl ude in this w o rk the photometry from 2009 
presented in iBarlow et all (|201Ch and a f ew light curves 
obtai ned with the Goodman Spectrograph (iClemens et al.l 
120041 1 on the 4.1-m SOAR telescope. Appendix [Al presents 
a detailed log of our photometric observations. 

2.1 Reduction &: Analysis 

Images were bias-subtracted and flat-fielded in IRAF0 
using standard procedures. We then extracted our pho- 
tometry using Antonio Kanaan's external IRAF pack- 
age ccdjisp. We chose aperture radii that maximized the 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation. 




Figure 1. Field image of CS 1246. Image shown is a stack of 
100 30-s Open-filter frames taken with PROMPT 3. CS 1246 
is marked with a star symbol, while the comparison stars used 
for the majority of the light curves are labelled C1-C3. The 
arrow indicates the drift direction and magnitude of the field 
as it appears in the frames over three hours of tracking with 
PROMPT 3. 

signal-to-noise ratio (S/N) in the nightly light curves and 
used sky annuli to estimate and subtract off counts from 
the sky. As the field is densely packed with stars, we were 
careful to avoid contamination of the apertures and an- 
nuli from neighboring stars. We removed sky transparency 
variations by dividing our light curves by those of con- 
stant comparison stars. Such stars were chosen by selecting 
from those in the field that were photometrically constant 
in time, not contaminated by nearby stars, substantially 
brighter than CS 1246, and not saturated during a 30-s 
exposure. Additionally, we avoided comparison stars that 
were likely to drift out of the frame during the PROMPT 
runs due to an approximately 0.4 arcmin hour -1 track- 
ing error in the mounts. These selection criteria provided 
us with three dependable comparison stars, as shown in 
the field image in Figure [T] Lastly, we fit and normalised 
the light curves with parabolas to remove residual atmo- 
spheric extinction effects. This process can remove real 
variations in the stellar flux on the order of the run lengths, 
making it difficult to detect lower-frequency j-mode os- 
cillations, if they exist. All time stamps were converted 
to a barycentric Juli an ephemeris date (BJED) u sing the 
WQED suite (v2.0; [Thompson fc Mullallvl [20091 ). which 
employs the method of IStumpfll (Il98pf h 

We analysed our light curves with a combination 
of Fourier analysis an d least-squares fits o f sine waves, 
using both Perio d04 (|Lenz fc Breeerl I2005T ) and MPFIT 
(|MarkwardtJl2009r ). Initially, the period, phase, and ampli- 
tude were left as free parameters in the fits, but once we 
obtained a better determination of the period (see Sj3j> , we 
refit all of the data with the period fixed. Figure [5] shows 
a representative light curve from one of our photometry 
runs above its amplitude spectrum. 

2.2 Night-by-night results 

Each of the individual 89 light curves is dominated by 
the 371. 7-s pulsati on mode (hereafter /i) reported by 
IBarlow et alfl|2010l l and shows no other significant signals 
in the p-mode range of the amplitude spectrum (1600- 
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Figure 2. (Top) Sample light curve of CS 1246 from 2010 Jan 
18. The data shown were obtained with PROMPT 3 over 2.9 
hours using 30-s exposures without a filter. (Bottom) Ampli- 
tude spectrum of the light curve, which is dominated by a 371.7- 
s signal. 
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Figure 3. (Top) Amplitudes of /i taken from the least-squares 
fits to the individual 2010 light curves. The data show a de- 
creasing trend of nearly 0.9 mma month , which is marked by 
a dotted line. (Bottom) Fourier transform of the nightly am- 
plitudes with the linear trend removed. The mean noise level 
is 0.25 mma. The dashed vertical line marks the location of a 
14.1-day oscillation, and the dotted horizontal line represents 
the amplitude at four times the mean noise level. 



14000 fiRz) after pre-whitening the signal. Low-amplitude 
signals appear in some of the nightly Fourier transforms 
(FTs) at frequencies near the g-mode regime, but we hes- 
itate to claim these signals as real since 1/f noise and im- 
properly removed extinction effects could easily account 
for those structures. 

Our least-squares fits show the amplitude of /i de- 
creased from 20 mm^flin 2010 Jan to 16 mma in 2010 May 
at a rate of nearly 0.9 mma month -1 . Figure [3] presents 
this trend above its FT. We exclude the 2009 data from 
this plot since they were taken primarily through different 
passbands from the 2010 data and the measured ampli- 
tude of A strongly depends on the observed wavelengths 
(|Barlow et al.| [2010). If the rate of decline remains linear, 
/i will have an amplitude below our detection limits (in 
nightly light curves) by 2011 Nov. As we discuss later, 
there is no detectable periodic signal in the amplitudes. 

The frequencies of /i in all the light curves agree 
within their error bars, which were typically around 1.2 
(iHz. To investigate smaller changes in the frequency as 
well as variations in the phase, we analysed our complete 
data set using the O-C technique, as discussed in the sec- 
tion that follows. 



3 THE O-C DIAGRAM 

We began our construction of the O-C diagranQ by calcu- 
lating a linear ephemeris for the times of maxima in the 
light curve (C) of the form 

C = T C + P C E, (1) 

where T c is a reference time of maximum, P c the pulsa- 
tion period, and E the cycle number measured from T c . 
Since keeping track of E correctly over an extended pe- 
riod of time requires an exceptionally accurate starting 
frequency, we used the period resulting from the least- 
squares fits to the combined 2010 light curve for P c . Our 
observing run on 2010 Mar 10 falls near the middle of this 
combined light curve, and, consequently, we employed the 
time of maximum determined from its least-squares fit for 
T c . The observed times of maxima (O) were taken from 
the least-squares fits to the individual light curves, and 
their corresponding cycle numbers were computed using 
equation |T|). 

Figure [4] presents the O-C diagram created by sub- 
tracting the calculated from the observed times of max- 
ima. A sinusoidal oscillation dominates the overall struc- 
ture. The sinusoid is superimposed on a parabolic trend 
indicative of a secular change in the pulsation period. The 
presence of the sinusoid is even more apparent in the FT 
of the O-C diagram with the quadratic term removed, as 



J Amplitudes are given in units of milli-modulation amplitude For a review of the basic principl es behind the O-C method, 

(mma), or parts-per-thousand. 10 mma corresponds to 1% we refer the reader to iKeplerl (1993). 
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Figure 4. (Top) The O-C diagram for CS 1246. O-C values were computed using f\ and a linear ephemeris. The diagram is 
dominated by a strong sinusoidal pattern with a period of 14.1 days overlaid on a parabola. (Middle) O-C points after removal of the 
quadratic term and (Bottom) after removal of both the parabola and sine wave. The mean noise level in the pre-whitened diagram 
is 0.75s. 



shown in Figure [S] To quantify these structures, we per- 
formed a simultaneous fit to the O-C values including both 
parabolic and sinusoidal terms using the expression 



O-C = AT - 



APE + ^PPE 2 



Aflin^ + tf). (2. 



We used the IDL routine MPFIT (|Markwardtl |2009h . 
which employs the Levenberg-Marquardt method, to per- 
form a non-linear least-squares fit of equation ([2]) to the 
data. During this process, the points were weighted by 
their phase error bars as determined from the least-squares 
fits. Inspecting both equations JTJ and ([2]), one finds that 
the resulting values for AT and AP provide corrections 
to our initial estimates for the reference time of maxi- 
mum and period from which we derive a final pulsational 
ephemeris for the times of maxima: 



= T + PE + \-PPE 2 + A sin 



( 2-kE 

v~fr 



(3) 



Table [T] displays the best-fit parameters for this ephemeris. 

The parabolic component of the fit indicates a sec- 
ular decrease in the pulsational period on the order of 1 
ms every 1.7 yrs. We attribute this variation to structural 
changes in CS 1246 as it evolves, as we discuss further 
in Sj6] Removing this signal from the O-C data (Figure 
U middle panel) more clearly reveals the two-week phase 
oscillation, which has a semi-amplitude of nearly lis. We 
point out that even the 2009 data, which are separated 
from the 2010 points by more than five months, phase 
well to this oscillation. We find no additional phase vari- 
ations after removing the parabolic and sinusoidal terms 



in either the O-C diagram (Figure 01 bottom panel) or its 
FT (Figure [5] bottom panel). The mean noise level in the 
FT of the pre-whitened O-C diagram is 0.75 s. 
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Table 1. Ephemeris parameters for the times of light maxima given in Eqn.[3] 



Parameter 


Value 


Error 


Units 


Comments 


T 


2455266.604605 


± 0.000004 


days 


time of light maximum 


P 


371.69162 


± 0.00004 


s 


fl pulsational period" 


P 


-1.9 x lO" 11 


± 0.3x 10- 11 


s s —1 


fl pulsational period change" 


A 


10.7 


± 0.4 


s 


fortnightly phase variation semi-amplitude 


11 


14.103 


± 0.010 


days 


fortnightly phase variation period 


4> 


0.023 


± 0.007 


cycles 


fortnightly phase variation phase" 



"as measured at T D 
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Figure 5. Fourier transform of the O-C values after removal of 
the quadratic term (middle panel). For reference, the window 
function for a sine wave sampled in the same manner as our 
data is also shown (top panel). The Fourier transform of the 
data with the parabola and sine wave removed (bottom panel) 
has a mean noise level of 0.75 s. The dashed line marks the 4-cr 
level. 



4 THE COMBINED LIGHT CURVE 

Three global variations associated with the pulsations of 
fi have been presented thus far: an amplitude decrease, a 
two- week phase oscillation, and a decrease in the pulsation 
period. Each of these phenomena should manifest itself in 
the FT of the combined light curve in a particular way. 
For this reason and to look for lower-amplitude signals 
not detectable in the nightly light curves, we combined 
the light curves and computed the Fourier transform of the 
result. The 2009 light curves were ignored in this analysis 
since they are few in number compared to the 2010 data 
and the large gaps they introduce to the combined light 
curve overly complicate the window function. 

Figure HJJi shows the amplitude spectrum of the com- 
plete light curve with an expanded view around fi , which 
clearly dominates the spectrum. After removing a sine 
wave of fixed frequency and amplitude from the light curve 
and re-computing the FT (Figure [6p), one can see signifi- 
cant power at the location of the first harmonic of the main 
mode (5381 fiHz) with an amplitude near 0.6 mma. The 
probability that random noise could produce this much 
signal by chance at exactly 2/i is less than 10 -12 ; conse- 



quently, we regard this signal as real. Evaluating the am- 
plitude spectra of light curves combined month-by-month 
provides further support for the presence of the first har- 
monic. 

Even more prominent than the first harmonic, the am- 
plitude spectrum of the pre-whitened light curve reveals a 
substantial amount of residual structure surrounding the 
location of fi . The most significant of these signals is an 
apparent 0.82-/j,Hz splitting of fi and its sidelobes, which 
is shown in the right half of Figure [6p. With an amplitude 
nearly 10% that of the original mode, these components 
are split by exactly 1/14.1 day, corresponding to the pe- 
riod of the phase oscillation observed in the O-C diagram. 
An additional splitting with an even tighter frequency sep- 
aration and much smaller amplitude is also visible in the 
panel. 

To test whether the above structures are consistent 
with the aforementioned three global variations associated 
with fi, we modeled the complete light curve with the 
function 



I(t) = (ai + d\t) sin 



+ffl2 sin 



2vrf 
P + Pt 
2-Kt 



- + -t 

L 2 2 



+ 4>i + <p(t) 

<j>2 + 2<p(t) 



(4) 



which includes a fundamental mode and first harmonic 
with linear period changes. Additionally, the fundamental 
mode is given a linear amplitude variation. Both terms 
contain a phase oscillation given by 



2-kA 
P + Pt' 



2-ut 



(5) 



Fixing P to its value determined from the O-C diagram, 
we let all other parameters vary freely and used MPFIT 
to find the best-fitting values shown in Table [2] We note 
that all fitted parameters determined previously in this 
manuscript (P, A, II, o'i, ffli, 02) using other methods (O- 
C diagram, night-by-night light curve fits) agree with those 
in the table to within the errors. Subtracting this fit from 
the data and re-calculating the FT (Figure [EJ:) shows our 
model removed essentially all of the signals centred around 
fi. The only signal remaining in the amplitude spectrum 
after this process is a low-amplitude peak offset by 15 /xHz 
from fi at 2715 /iHz. The frequency separation from fi 
implies an interaction timescale of 0.5 day. The origin of 
this signal is unclear, but the probability that noise could 
produce it by chance is less than 10 -5 . 

Other scenarios could also lead to the residual FT 
structure observed around fi in Figure |Bp, but they would 
have to do by conspiracy what our proposed model does 
naturally. The equidistant triplets might be produced by 
(a) the splitting of fi by stellar rotation, (b) the presence 
of additional independent pulsation frequencies symmet- 
rically disposed about the original, (c) one independent 
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Figure 6. (a) Amplitude spectrum of the combined 2010 light curve (left column) with expansions around fx (right column), (b) 
Spectrum of the light curve pre- whitened by a sine wave with frequency fx, which reveals a 0.82-/iHz splitting of fx. Note also the 
presence of the first harmonic to fx ■ (c) Amplitude spectrum of the data pre- whitened by the best fit of equation Q . The mean 
noise level is around 0.074 mma. 



pulsation frequency and a combination mode at the differ- 
ence between it and the first harmonic, and (d) amplitude 
variability of fx. None of these is appealing nor as effec- 
tive as the model of equation (U) at removing the residual 
power. In particular, all of the alternative models predict 
an amplitude modulation of fx with a period of two weeks, 
which is not seen in the Fourier transform of the ampli- 
tude fits after the decreasing amplitude trend is removed 
(Figure [3] bottom panel). 



5 BINARY SYSTEM PARAMETERS 

Under the assumption that the O-C variations represent 
reflex motion due to an orbital companion, we have fitted 
the best orbital parameters. Using the fits in Table [1] CS 
1246 moves about a barycentre at least 11 light-seconds 
(4.7 R.0) away with a period of two weeks. If the orbit is 
perfectly circular, the phase variation will be sinusoidal; 
otherwise, the shape will deviate from a sine wave. As 
some s ubdwarf binaries show small eccentricities in their 
orbits (|Edelmann et al.l l2006). we decided to evaluate the 
eccentricity of the orbit by investigating both the phase- 
folded O-C curve and the Fourier transform thereof. Figure 
presents the phase-folded O-C diagram after removal of 
the parabolic trend. 

The best-fitting orbital parameters to the folded O-C 
diagram yield an eccentricity of e = 0.04 ± 0.04, which we 
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Figure 7. (Top) Phase-folded O-C diagram with parabola re- 
moved. The O-C points were folded on the fortnightly peri- 
odicity and averaged into ten bins. The dotted line shows the 
circular orbit fit to the O-C diagram. (Bottom) Residual O-C 
values after removal of a sinusoid fit. 



take to be consistent with a circular orbit. The absence of 
a prominent first harmonic (at 0.14 day" 1 ) in the FT of 
the O-C diagram (Figure [S] bottom panel), which would 
be expected for a considerably elliptical orbit, corrobo- 
rates this conclusion. Some degree of eccentricity might 
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Table 2. Best fit parameters to Eqn. {1} 



Parameter 


Value 


Error 


Units 


Comment 


"i 


17.71 


± 0.075 


mma 


fl amplitude" 


a'l 


-0.0298 


± 0.0027 


mma day 


fl amplitude change 6 


a 2 


0.5858 


± 0.075 


mma 


2/i amplitude a,c 


P 


371.691646 


± 0.000033 


s 


fl period" 


P 


-1.88 


(fixed) 


10 -11 s s" 1 


fl period changc d 


A 


10.69 


± 0.35 


s 


fl phase oscillation semi-amplitude 


11 


14.067 


± 0.032 


day 


fl phase oscillation period 


01 


0.85654 


± 0.00070 


cycles 


fl phase" 


<p2 


0.435 


± 0.020 


cycles 


fl phase" 




0.8456 


± 0.0055 


cycles 


phase of the fi phase oscillation" 



" at t=0 corresponding to BJED 2455271.197436 

b from 2010 Jan to 2010 May 

c with f2 fixed to 2fi 

d fixed to value from Tabled 



Table 3. System Parameters 



Param Value Error Units Comment 



n 14.105 ± 0.011 days orbital period 

K 16.6 dt 0.6 km s 1 RV scmi-amplitude a 

/ 0.0066 ± 0.0007 Mg mass function 

a 0.0910 ± 0.0003 AU separation distance 6 

0.0963 ± 0.0003 AU separation distance c 

m sin i 0.115 dt 0.005 ^0 minimum companion mass ' 

0.129 i 0.005 ^0 minimum companion mass c ' 



an eclipse only rules out inclination angles greater than 
89 deg, which increases the minimum mass we report by a 
negligible amount. 

Under the assumption of randomly distributed orbital 
plane orientations, one computes the probability of observ- 
ing a system at an inclination angle i less than or equal to 
6 as 

Probability (i < 9) = 1 - cos 9 (7) 

From this we calculate a 96% probability that the compan- 
ion mass is less than 0.45 Mq, hence the object is most 
likely a low-mass white dwarf or late- type main sequence 
M-dwarf. Optical spectra of CS 1246 (|Barlow et aljhoid, 
show no clear signatures of a companion, but this result 
is consistent with both companion possibilities. Optical 
reflection effects are often observed in sdB+dM systems, 
but even if the albedo of the M-dwarf were 1.0, the rel- 
atively large separation distance would result in a maxi- 
mum flux increase less than 0.01 mma, a modulation well 
below our detection limits. An infrared excess would be 
expected for CS 1246 if the companio n is a main sequenc e 
star with mass exceeding 0.45 Mq (|Lisker et al.l [2005h - 
Unfortunately, the presence of the Coalsack Nebula along 
our line of sight significantly reddens the system flux and 
complicates identification of an infrared excess from a com- 
panion. 

6 THE EVOLUTION OF CS 1246 

Clues to the current evolutionary state of CS 1246 may be 
found through cons ideration of the pulsati onal P. Accord- 
ing to the models of lCharpinet et all l|2002h , changes in the 
periods of acoustic modes in EHB stars are dominated by 
two primary phases. During the first, a decrease in the sur- 
face gravity causes the periods to increase steadily as the 
sdB star evolves away from the zero-age extended horizon- 
tal branch (ZAEHB). Model stars display positive values 
of P in this phase. The second phase of evolution begins 
approximately 90 Myr after the ZAEHB when the den- 
sity of thermonuclear fuel in the core begins to decrease 
dramatically. Responding to the reduced energy produc- 
tion, the model stars contract, and the resulting increase 
in surface gravity leads to a decrease in acoustic mode pe- 
riods and hence a negative value for P. This phase lasts 
until the onset of the post-EHB phase approximately 110 
Myr after the ZAEHB, at which point the He in the core 
is completely exhausted. Thereafter, He fusion may be- 



a circular orbit approximation 

assuming an sdB mass of 0.39 Mq 
c assuming the canonical sdB mass of 0.47 Mq 
rf assumes no error bar on the sdB mass 

be expected since the circularization timescale for a bi- 
nary with a two week orb ital period exceeds 2 billion years 
(|Tassoul fc Tassoulll 19921 ) but would not significantly alter 
the results of our mass determination. 

Without information concerning the orbital motion 
of the unseen companion, we are limited to computing the 
mass function for the system: 

_ m 3 sin 3 i _ PK 3 
' ~ (m + M) 2 ~ 2ttG () 

The companion mass (m), sdB mass (M), and inclination 
angle (i) are free parameters, while the period (P) and 
RV semi-amplitude (K) of the sdB are taken from our fits 
to the O-C diagram. To compute the mas s of the com- 
panio n, we must assume a mass for the sdB. lBarlow et all 
(2010) showed that under the assumption of a radial oscil- 
lation, a mass of 0.39 to ls M could be derived using the 
Baade-Wesselink method. However, since the error bars 
are large and the radial nature of fi remains unproven, 
we u se both this mass and the canonical sdB mass (0.47 
M0: lHan et al.li2002l . l2l)03h to derive the system parame- 
ters shown in Table 

The minimum mass derived for the unseen compan- 
ion is 0.12 Mq. Further limitations could be placed on this 
mass by constraining the inclination angle of the system. 
If the orbital plane is nearly edge-on (i ~ 90 deg), primary 
and secondary eclipses will occur when the O-C variations 
are close to their maximum and minimum values, respec- 
tively. To look for an eclipse, we carefully monitored the 
light curve using three of the PROMPT telescopes on 2010 
Jun 21 (UTC) at a time when equation ([3| predicts a pri- 
mary eclipse of the sdB in the case of an edge-on orbital 
plane. None of the light curves showed any sign of a signal 
loss. Given the large separation distance, the absence of 
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gin in the shell, the radius begins to increase, and the 
star heads towards the white dwarf cooling sequence after 
passing through the sdO regime of the H-R diagram. Most 
helium-deficient sdO stars are believed to be the progeny 
of sdB stars evolving in this way. 

Our negative value of P, interpreted as an evolution- 
ary effect, implies CS 1246 is nearing the end of its life 
as an EHB star or has already done so. The contraction 
resulting from the depletion of He in the core yields an 
increased surface gravity, which leads to a decrease in the 
pulsat ional period of -571 ± 85 s Mvr" 1 . fcharpinet et al.l 
(2002) calculated values of P for sdB acoustic modes (their 
Appendix C) for representative models at different ages. 
Comparing our P to all of their theoretical ones, we find 
our rate of change to be more rapid than all of their listed 
values but most closely matches that of a radial mode in 
a model with an age of 106.58 Myr and a hydrogen layer 
mass of 0.0042 Mq. Our faster P could result from CS 
1246 having a thicker hydrogen lay er than any of their 
models, since ICharpinet et al.l (|2002r i show increasing the 
hydrogen layer mass leads to faster period changes. It is 
also possible CS 1246 exists in an evolutiona r y epo ch not 
represented in their table. ICharpinet et alJ l|2002t ) note 
that in advanced stages of EHB evolution, when the pe- 
riod decreases most rapidly, they used simple forward dif- 
ferencing to estimate theoretical P values instead of cubic 
spline fittings, which they used in earlier EHB stages when 
the periods are more stable. A potential consequence of 
this differencing is the over- or under-estimation of period 
changes at these advanced evolutionary stages. 

Having stated these possibilities, we stress that our 
measurement of P is at most an upper limit to the evo- 
lutionary rate of CS 1246. Secular processes in sdB stars 
work on slower timescales than other processes affecting 
the phase or frequency of a mode, whic h in some case s 
drive the observed P to faster changes l|Kawalerl l2010h . 
Proper motion, for example, can contrib ute to the mea - 
sured period changes in pulsating stars (|PaidosdlT995r ). 
but the expected contribution in our case is several orders 
of magnitude smaller than our observed P and inconse- 
quential to our analysis. 



7 DISCUSSION & OUTLOOK 

Using the O-C technique, we have discovered an apparent 
orbital reflex motion in CS 1246 and measured an upper 
limit of P for its main pulsation mode. This marks the sec- 
ond time such quantities have been derived for an sdBV 
star using thi s meth od. The first measurements were from 
ISilvotti et al.l (2007). The measured rate of period change, 
if secular, implies the star has either already exhausted 
the He in its core or is quickly approaching the post-EHB 
stage. Models show the structure of sdB stars in this evo- 
lutionary state change dramatically fast, and it may even 
be possible to measure P within a few years. 

An analysis using Kepler's law shows the companion 
has a minimum mass near 0.12 Mq and is most likely a 
white dwarf or late-type main sequence star. The separa- 
tion distance implies the companion was potentially inside 
the envelope of the sdB progenitor during its red giant 
phase and might even be responsible f or the ejection of 
the envelope fsee lHan et al.l l2002. 200^)- We cannot rule 
out the presence of additional bodies in the system, but we 
can limit the possible combinations of companion masses 
and separation distances using the detection limits from 
the Fourier transform of our O-C data. The shaded region 
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Figure 8. Approximate sensitivity of our O-C diagram to com- 
panions around CS 1246. The shaded region marks objects that 
are not detectable using the O-C data shown in FigurefJ] Com- 
panions falling within the white region could be detected in our 
data set, with certain restrictions on their orbital inclination an- 
gle. For these calculations, we assume a mass of 0.47 Mq for CS 
1246. The two dashed lines mark the combinations of compan- 
ion masses and separation distances leading to orbital periods 
equal to the longest (right dashed line) and shortest (left dashed 
line) periods detectable using our O-C data, as computed us- 
ing Kepler's third law. The solid line shows the combinations 
of masses and distances resulting in reflex motions equal to our 
detection limits in the Fourier transform of the O-C diagram, 
under the assumption of an edge-on orbital plane. Masses be- 
low this line induce reflex motion below our detection limits, 
regardless of their inclination angles. The star symbol marks 
the location of the companion we detected (plotted using the 
minimum mass), and the two dotted horizontal lines show the 
boundaries between planet and brown dwarf masses (lower dot- 
ted line) and between brown dwarf and stellar masses (upper 
dotted line). 



in Figure [8] shows combinations of masses and separation 
distances leading to phase wobbles that are undetectable 
using our current data set. Our ability to detect a com- 
panion with orbital parameters inside the white region 
depends on the inclination angle of its orbit. The lower 
solid line shows the combinations of system parameters 
resulting in reflex motions equal to our detection limits, 
under the assumption of an edge-on orbital plane; the line 
moves up with increasing inclination angle. The region is 
bounded on the left and right by dashed lines correspond- 
ing to the minimum and maximum orbital periods we can 
detect given the sampling rate and baseline of our O-C 
data; thus, we could expand the area of the white region 
by obtaining phase measurements more frequently or by 
monitoring the system for a longer period of time. Regard- 
less of the sampling, additional phase measurements will 
lower our detection limits and move the lower boundary 
(solid line) to smaller masses. As a reference, the star sym- 
bol shows the location of the single companion we detect 
around CS 1246 (plotted with its minimum mass). If the 
observed phase oscillation is the result of a binary com- 
panion, CS 1246 would show the second-lowest RV semi- 
amplitude and the third-longest period of the ~85 sdB 
binary syst ems for which orb ital parameters are known 
(Table A.l. lGeier et"afll2011cft . 

More statistics are needed to draw conclusions about 
the formation channels leading to sdB stars from the or- 
bital parameters of binary systems. Searches for radial ve- 
locity variations, as are being carried out by the MUCH- 
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FUSS survey l|Geier et all 1201 id ), should uncover a wide 
variety of sdB binary systems, although they are some- 
what biased towards higher-mass systems with shorter or- 
bital periods. Pulsating sdB stars may also be used to 
find unseen compani ons using the O-C di agram, as shown 
in this work a nd by ISilvotti et all (120071 ). and the EXO- 
TIME project (|Schuh et alJl201Ch is currently building up 
O-C diagrams for several sdBV stars to look for oscilla- 
tions in the pulse timings. Although both of these surveys 
will be effective at uncovering binary companions to sdB 
stars, their sampling rates differ significantly from the one 
presented in this study, and, consequently, they will not 
be sensitive to the same parameter space shown in Figure 

El 

Our result provides an important link between the two 
aforementioned methods for detecting unseen companions 
since an observable radial velocity variation should accom- 
pany our fortnightly phase oscillation. The O-C diagram 
is a significant tool used to discover and characterize the 
nature of unseen companions, including extrasolar planets, 
but no corroborative RV variations have been measured to 
date. The O-C oscillation of CS 1246 implies a two-week 
velocity variation with a semi-amplitude well within the 
range of medium-to-high resolution spectrographs. Conse- 
quently, we are planning follow-up observations to mea- 
sure the RV variations. If the expected RV variation is ob- 
served, it will be the first independent confirmation that 
it is possible to detect stellar and planetary companions 
using the O-C technique and sdB pulsations as the clock. 
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Table Al. Observations Log 



Date 


Start Time 




^ cycl e 


Length 


Filter 


Scope" 


Date 




Start Time 




T cycle 
(s) 


Length 


Filter 


Scope 


(UTCJ 


(UTC) 


0) 


0>) 


(hr) 






(UTC) 


(UTC) 


(s) 


(hr) 






2009-03-29 


05:13:09 


40 


46 


1.23 


V 


P4 


2010-02 


21 


02:00:07 


30 


36 


2.91 


Open 


P3 


2009-04-01 


00:51:02 


40 


46 


1.40 


V 


P5 


2010-02 


22 


02:07:14 


30 


36 


2.41 


Open 


P3 


2009-04-03 


01:16:16 


20 


26 


1.17 


Open 


P4 


2010-02 


23 


01:52:20 


31 


37 


2.48 


Open 


P3 


2009-04-03 


23:59:03 


80 


86 


9.43 


u' 


P3 


2010-02 


25 


02:27:30 


30 


36 


2.40 


Open 


P3 


2009-04-04 


00:57:03 


40 


52 


9.57 


g' 


P2 


2010-02 


27 


01:37:09 


30 


36 


3.11 


Open 


P3 


2009-04-04 


00:01:48 


40 


47 


9.57 


r ' 


P4 


2010-02 


28 


08:04:58 


20 


26 


1. 78 


S8612 


SOAR 


2009-04-04 


00:01:52 


40 


47 


8.66 


1 


P5 


2010-03 


01 


02:37: 10 


20 


26 


6.79 


B 


SOAR. 


2009-04-04 


23:57:56 


80 


86 


8.61 


u' 


P3 


2010-03 


05 


01: 13:28 


30 


36 


2.41 


Open 


P3 


2009-04-04 


23:57:06 


40 


52 


9.31 


g' 


P2 


2010-03 


06 


01:09:45 


30 


36 


1.61 


Open 


P3 


2009-04-04 


23:57:42 


40 


47 


9.29 


y 


P4 


2010-03 


07 


01:06:47 


30 


36 


2.92 


Open 


P3 


2009-04-04 


23:57:39 


40 


47 


9.43 


1 


Po 


2010-03 


11 


00:44:23 


30 


36 


3.20 


Open 


P3 


2009-04-05 


23:56:42 


80 


86 


9.70 


u' 


P3 


2010-03 


13 


02:09:46 


30 


36 


1 .64 


Open 


P3 


2009-04-05 


23:55:47 


40 


52 


9.69 


g' 


P2 


2010-03 


14 


00:32:23 


30 


36 


3. 19 


Open 


P3 


2009-04-05 


23:56:36 


40 


47 


9.68 


r ' 


P4 


2010-03 


16 


01:00:50 


30 


36 


2.59 


Open 


P3 


2009-04-05 


23:56:32 


40 


47 


9.69 


*' 


P5 


2010-03 


18 


01:21:58 


30 


36 


4.37 


Open 


P3 


2009-04-16 


09:07:06 


10 


16 


1.10 


R 


SOAR 


2010-03 


21 


06:02: 18 


30 


36 


3.41 


Open 


P3 


2009-04-17 


01:00:22 


80 


86 


6.14 


350-620nm b 


SOAR 


2010-03 


23 


01:20:24 


30 


36 


2.56 


Open 


P3 


2009-05-14 


05:09:12 


20 


26 


3.08 


V 


P5 


2010-03 


27 


03:00:52 


30 


36 


2.14 


Open 


P3 


2009-05-15 


00:57:56 


20 


26 


5.45 


V 


P5 


2010-03 


28 


02:51:33 


30 


36 


1.71 


Open 


P3 


2009-07-23 


23:36:05 


40 


46 


3.25 


V 


P4 


2010-03 


29 


06:30:50 


30 


36 


2.99 


Open 


P3 


2009-08-02 


23:35:20 


40 


46 


2.67 


r prime 


P3 


2010-03 


30 


00:28:42 


30 


36 


3.00 


Open 


P3 


2009-08-04 


00:00:41 


40 


46 


1.48 


rprime 


P3 


2010-03 


31 


05:06:29 


30 


36 


2.99 


Open 


P3 


2009-08-05 


00:00:48 


40 


46 


2.12 


r prime 


P3 


2010-04 


01 


00:29:16 


30 


36 


1.98 


Open 


P3 


2010-01-17 


04:18:07 


30 


36 


3.25 


Open 


P3 


2010-04 


02 


00:17:15 


30 


36 


3.97 


Open 


P3 


2010-01-18 


04:13:15 


30 


36 


2.90 


Open 


P3 


2010-04 


03 


00:50:48 


30 


36 


3.49 


Open 


P3 


2010-01-19 


04:09:46 


30 


36 


3.88 


Open 


P3 


2010-04 


04 


23:58:13 


30 


36 


3.08 


Open 


P3 


2010-01-20 


04:05:29 


30 


36 


3.94 


Open 


P3 


2010-04 


06 


01:26:05 


30 


36 


7.95 


Open 


P3 


2010-01-21 


04:01:50 


30 


36 


3.90 


Open 


P3 


2010-04 


07 


00:46:28 


30 


36 


5.28 


Open 


P3 


2010-01-22 


04:52:00 


33 


39 


2.91 


Open 


P3 


2010-04 


08 


00:41:50 


30 


36 


4.67 


Open 


P3 


2010-01-23 


03:53:38 


25 


31 


3.38 


Open 


P3 


2010-04 


09 


00:42:05 


30 


36 


2.98 


Open 


P3 


2010-01-24 


04:39:12 


30 


36 


3.42 


Open 


P3 


2010-04 


10 


00:38:21 


30 


36 


3.59 


Open 


P3 


2010-01-26 


03:42:18 


30 


36 


1.48 


Open 


P3 


2010-04 


12 


00:12:12 


30 


36 


5.56 


Open 


P3 


2010-01-29 


01:51:06 


30 


36 


3.74 


Open 


P3 


2010-04 




00:07:47 


30 


36 


4.99 


Open 


P3 


2010-01-31 


03:22:26 


30 


36 


3.93 


Open 


P3 


2010-04 


14 


00:03:33 


30 


36 


3.70 


Open 


P3 


2010-02-01 


03:18:25 


30 


36 


4.96 


Open 


P3 


2010-04 


15 


00:03:48 


30 


36 


1.90 


Open 


P3 


2010-02-02 


06:41:19 


30 


36 


1.93 


Open 


P3 


2010-04 


19 


03:51:05 


30 


36 


4.29 


Open 


P3 


2010-02-04 


03:06:53 


30 


36 


4.59 


Open 


P3 


2010-04 


19 


23:49:39 


30 


36 


5.41 


Open 


P3 


2010-02-07 


01:56:24 


30 


36 


2.69 


Open 


P3 


2010-04 


20 


23:39:39 


30 


36 


7.40 


Open 


P3 


2010-02-08 


04:19:12 


30 


36 


5.38 


Open 


P3 


2010-04 


22 


23:37:12 


30 


36 


6.63 


Open 


P3 


2010-02-09 


04:00:47 


30 


36 


3.84 


Open 


P3 


2010-04 


24 


03:47:37 


30 


36 


1.84 


Open 


P3 


2010-02-10 


04:01:02 


30 


36 


4.04 


Open 


P3 


2010-04 


25 


22:58:30 


30 


36 


2.82 


Open 


P3 


2010-02-13 


03:49:45 


30 


36 


2.55 


Open 


P3 


2010-04 


27 


03:31:02 


30 


36 


3.66 


Open 


P3 


2010-02-14 


03:59:15 


30 


36 


2.33 


Open 


P3 


2010-04 


27 


23:33:40 


30 


36 


2.50 


Open 


P3 


2010-02-16 


03:51:20 


30 


36 


2.33 


Open 


P3 


2010-04 


30 


23:16:38 


30 


36 


3.58 


Open 


P3 


2010-02-17 


02:49:32 


30 


36 


3.00 


Open 


P3 


2010-05 


03 


02:56:35 


30 


36 


4.16 


Open 


P3 


2010-02-18 


02:11:52 


30 


36 


3.11 


Open 


P3 


2010-05 


04 


23:28:40 


30 


36 


3.11 


Open 


P3 


2010-02-19 


02:08:17 


30 


36 


3.23 


Open 


P3 


2010-05 


07 


04:31:20 


30 


36 


2.92 


Open 


P3 


2010-02-20 


02:58:25 


30 
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a u p „ _ PROMPT 

^light curve produced from time-series spectra over the wavelength range 3500-6200 A 
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The O-C diagram of CS 1246 



Table Bl. Observed times of light maxima and corresponding O-C values 
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